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I  INTRODUCTION 


The  mechanical  properties  of  metals  have  long  been  of  interest  because  of  their  use  as 
structural  materials.  Even  so,  much  of  the  understanding  of  the  mechanical  properties  of 
metals  is  from  empirical  data  and  macroscopic  models.  In  semiconductors,  structural  proper¬ 
ties  are  also  of  interest  because  of  their  relationship  to  the  material’s  manufacturability,  as  well 
as  the  impact  that  structural  defects,  such  as  dislocations  and  vacancies,  may  have  on  electronic 
properties.  Thus,  we  strive  for  a  microscopic  understanding  of  the  mechanical/structural  pro¬ 
perties  of  the  semiconductors,  by  relating  the  observed  macroscopic  structural  properties  to 
microscopic  quantities  such  as  bond  lengths  and  energies,  elastic  constants,  and  so  on.  It  is 
only  through  a  microscopic  understanding  of  the  mechanical  behavior  of  semiconductors  that 
we  can  fully  understand  the  driving  mechanisms  of  structurally  inferior  materials  and  devise 
new  ways  to  improve  their  characteristics — for  example,  by  alloying,  varying  the  growth  tech¬ 
niques,  and  so  on.  Additionally,  new  nonequilibrium  growth  techniques  of  semiconductor 
alloys  may  result  in  material  possessing  long-  or  short-range  ordering,  which  in  turn  will 
influence  the  mechanical  properties  of  these  alloys.  Only  from  a  microscopic  theory  of  the 
*  mechanical  properties  of  the  materials  can  the  mechanical  properties  of  these  novel  new  materi¬ 

als  be  understood. 

The  remainder  of  this  report  is  organized  as  follows:  In  Section  II,  our  work  relating  to 
correlations  and  ordering  in  semiconductor  alloys  is  summarized.  In  Section  III,  results  of  the 
calculations  of  vacancy  formation  energies  and  surface  segregation  in  zinc-blende  semi¬ 
conductor  compounds  and  alloys  are  summarized.  Our  woric  on  the  elastic  constants,  hardness, 
and  cleavage  is  summarized  in  Section  IV.  Note  that  only  a  brief  summary  of  the  accom¬ 
plished  work  is  presented  in  cases  where  details  can  be  found  in  our  published  materials, 
which  are  listed  as  references  following  the  conclusion  of  the  text.  A  bibliography  of  publica- 
tioas  supported  in  whole  or  in  part  by  this  contract  follows  the  references. 
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II  CORRELATIONS  AND  ORDERING  IN  SEMICONDUCTOR  ALLOYS 


The  state  of  order  or  disorder  in  semiconductor  alloys  and  its  dependence  on  growth 
method  and  growth  and  processing  parameters  is  important  in  understanding  the  structural  pro¬ 
perties  of  these  materials.  The  first  step  towards  such  an  understanding  is  therefore  a 
knowledge  of  the  driving  forces  for  order  and  disorder.  We  have  developed  a  quasichemical 
approximation  to  study  effects  that  drive  short-range  order  in  the  pseudobinary  alloys.  The 
complementary  situation  of  alloys  with  long-range  order  has  also  been  calculated.  The  calcula¬ 
tions  for  the  short-  and  long-range  ordered  alloys  include  contributions  to  the  total  energy  from 
chemical  and  strain  terms.  One  further  contribution  to  the  total  energy  comes  from  Coulomb 
energies  in  the  alloys,  which  may  differ  substantially  in  the  ordered  structures  and  the  random 
alloys.  Our  calculations  show  that  the  Coulomb  energy  contribution  to  the  mixing  energies  is, 
in  general,  small,  although  in  some  lattice  matched  systems  where  the  usual  dominant  strain 
energy  is  small.  Coulomb  energies  may  dominate.  What  follows  is  a  brief  summary  of  our 
work  in  the  area.  General  conclusions  are  included  in  Section  V. 

A.  Quasichemical  Approximation  in  Pseudobinary  Alloys 

The  quasichemical  approximation  for  an  Ai_xBx  alloy  has  been  formulated  for  an  arbitrary 
lattice  and  choice  of  cluster.1  The  statistical  problem  of  the  average  number  of  each  class  of 
clusters  is  collapsed  into  a  polynomial  equation.  An  n-atom  cluster  of  type  j  is  characterized 
by  an  excess  energy  e-,  the  number  of  B  atoms  nj(B),  and  the  degeneracy  gr  If  Ej  is  a  linear 
function  of  n^B),  and  gj  is  a  binomial  coefficient  of  n  and  n^B),  then  the  cluster  populations 
will  be  random.  Strains  due  to  lattice-size  mismatches,  chemical  (electron-ion  interaction) 
differences,  and  differences  between  the  electron-electron  Coulomb  interactions  of  the  alloy 
constituents  will  drive  nonlinear  variations  of  Ej  on  nj(B).  The  gj  is  modified  by  coherent,  eter¬ 
nally  applied  stresses  and  temperature  gradients  present  during  crystal  growth.  The  conditions 
under  which  compounds  are  formed  or  spinodal  decomposition  occurs  have  been  derived.  The 
possibility  of  materials  consisting  of  arrays  of  two  kinds  of  domains,  one  a  random  alloy  and 
the  other  an  ordered  compound  have  been  considered.  The  theory  has  been  specialized  to 
pseudobinary  semiconductor  alloys  A]_XBXC  in  a  distorted  zinc-blende  structure.  Results  for 
many  systems  have  been  published.1-6 


The  generalized  quasichemical  approximation  we  have  developed  has  been  used  to  study 
the  phase  diagram  of  several  pseudobinary  alloys.5,7  The  solid  mixing  free  energies  were  cal¬ 
culated  from  the  generalized  quasichemical  approximation,1  while  the  liquid  mixing  free- 
energies  were  constructed  using  the  quasichemical  approximation  method  of  Stringfellow  and 
Greene.8  Knowing  the  free-energies,  the  liquidus  and  solidus  concentrations  for  a  given  tem¬ 
perature  were  calculated  by  the  common  tangent  method  equivalent  to  matching  the  chemical 
potential  for  the  liquid  and  solid  phases.  The  calculated  liquidus  and  solidus  curves  for  the 
alloys  considered,  GalnP,  GalnAs,  InSbAs,  HgCdTe,  HgZnTe,  and  CdZnTe  are  in  good  agree¬ 
ment  with  experiment. 

B.  Binding  Energies  and  Bond  Lengths  in  Ordered  Structures 

At  least  four  pseudobinary  alloys,  AlGaAs,  GalnAs,  GalnP,  and  GaAsSb  grown  from 
MBE  and  MOCVD  have  been  found  to  have  the  long-range  ordering  (LRO).  From  a  technol¬ 
ogy  point  of  view,  the  LRO  has  not  produced  unusually  electrical  or  mechanical  properties  that 
deserve  special  attention,  mainly  because  the  degree  of  LRO  is  rather  weak  and  the  ordering 
only  occurs  under  very  special  growth  conditions.  However,  the  LRO  may  provide  some  clue 
for  a  better  understanding  of  growth  processes  and  the  structural  properties  of  these  epitaxially 
grown  systems  as  compared  to  the  bulk  alloys. 

Whether  or  not  the  alloys  with  LRO  are  in  thermal  equilibrium  states  may  be  difficult  to 
test.  One  possible  test  is  a  simple  annealing  experiment.  If  the  LRO  is  restored  after  heating 
to  temperatures  above  the  transition  temperature  and  then  cooled  down  slowly,  the  alloy  will 
be  in  an  equilibrium  state.  If  not,  the  LRO  may  be  metastable.  The  difficulty  in  this  experi¬ 
ment  is  the  rates  of  diffusion  at  the  temperature  near  the  order-disorder  transition.  If  diffusion 
rates  are  slow  here,  the  time  necessary  to  restore  the  LRO  may  be  extremely  long  (effectively 
infinite).  An  accurate  calculation  of  the  total  energies  can  give  an  insight  into  this  problem. 
Below,  we  describe  our  unpublished  results  to  date  for  such  a  calculation. 

Our  calculation  is  based  on  the  first-neighbor  tight-binding  model.  The  total  energy  is  the 
sum  of  the  electron  energies  in  the  valence  bands,  Ee|,  and  the  nearest-neighbor  repulsive  pair 
energies,  V0.  There  are  two  major  differences  between  our  calculation  and  Harrison’s  bond- 
orbital  model.  First,  we  have  to  calculate  the  band  structures  for  the  appropriate  structures, 
i.e.,  zinc  blende  (zb),  copper-gold  I  (CuAu  I),  and  the  chacopyrite  (cp)  structures,  and  then  per¬ 
form  the  Brillouin-zone  summation  (BZ)  to  obtain  Ee).  The  second  difference  is  that  we  have 
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adjusted  the  values  and  the  scaling  rules  for  the  interaction  parameters  and  V0  to  improve  the 
agreement  between  the  calculated  quantities  and  the  experimental  values  for  the  constituent 
compounds  as  listed  in  Table  1 . 


With  four  s  and  p  orbitals  per  atom,  the  sizes  of  the  Hamiltonian  matrix  for  a  given  k  in 
the  BZ  are,  16  x  16,  and  32  x  32  respectively  for  the  CuAu  I  and  cp  structures.  A  uniform 
grid  of  125  K  points  inside  the  BZ  are  used  to  obtain  Ee].  The  total  energy  for  a  given  system 
is  calculated  at  four  guessed  lattice  parameters  a,  or  a  given  a  at  four  internal  distortion 
parameters;  these  results  are  used  to  locate  the  equilibrium  bond  lengths,  d!  and  d2,  and  the 
minimum  energy.  As  an  example,  Figure  1  shows  the  energy  per  bond,  1^,  as  a  function  of 
the  average  bond  length,  d,  for  GaAs,  InAs,  and  InGaAs  calculated  in  the  CuAu  I  structure. 
Table  2  lists  the  equilibrium  values  of  d,  dj  d2i  Eb,  two  kinds  of  excess  energies,  AE  and  AE\ 
to  be  discussed  below,  and  the  bulk  moduli  B  for  three  ordered  ternary  alloys  in  the  CuAu  I 
and  cp  structures.  Some  key  results  are  summarized  as  follows: 

(1)  The  average  d  =  (d!  +  d^/2  obeys  Vegard’s  law  very  well. 

(2)  The  calculated  d,  and  d2  fall  in  the  range  of  experimental  bimodal  distributions  found 
for  the  bulk  pseudobinary  alloys. 

(3)  The  alloy  bulk  moduli  are  a  few  percent  smaller  than  the  average  values  of  the  con¬ 
stituent  compounds. 

(4)  All  the  ternary  alloys  have  positive  excess  energies,  defined  as  AE  =  4E,, 
(ABC)  -  2Eb  (BC),  where  Eb  (AC)  and  Eb  (BC)  are  the  bond  energies  of  the  consti¬ 
tuent  compounds.  The  AE  values  are  larger  in  the  CuAu  I  than  in  the  cp  structures 
for  the  cation-substituted  alloys  such  as  InGaAs,  the  opposite  is  found  for  GaAsSb. 

(5)  All  ordered  alloys  have  very  substantial  negative  excess  energies,  AE',  which  is 
defined  similarly  to  the  AE  except  that  energies  of  the  two  constituent  compounds  are 
evaluated  at  a  lattice  constant  constrained  to  be  that  of  the  substrates. 

Table  1 

CALCULATED  AND  EXPERIMENTAL*  VALUES 
FOR  SIX  REPRESENTATIVE  COMPOUNDS 


AlAs 

GaP 

GaAs 

GaSb 

InP 

InAs 

d  (A) 

2.453 

2.360 

2.450 

2.642 

2.541 

2.625 

Eb  (eV) 

-1.8840 

(-L89) 

-1.7722 

(-1.78) 

-1.6239 

(-1.63) 

-1.4764 

(-1.48) 

-1.7406 

(-1.74) 

-1.5421 

(-1.55) 

B  (1011  dynes/cm2) 

mm 

8.836 

(8.87) 

7.555 

(7.48) 

5.919 

(5.63) 

1 

1 

5.823 

(5.80) 

♦Experimental  values  are  shown  in  parenthesis. 
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FIGURE  1  BOND  LENGTH  DEPENDENCE  OF  BOND  ENERGY 
IN  GaAs,  InAs,  AND  ln0.5Ga0  5As  WITH  THE  CuAu  I 
STRUCTURE 


The  last  two  results  are  the  most  relevant  to  the  order-disorder  transitions.  For  a  bulk 
alloy,  the  AE  is  a  measure  of  the  mixing  enthalpy,  and  the  disordered  alloy  will  be  immiscible 
below  the  critical  temperature,  Tc  =  AE/k.  Our  results  are  in  agreement  with  what  we  found  in 
our  detailed  studies  of  the  statistics  and  phase-diagrams  of  the  bulk  alloys.  However,  for  an 
epitaxially  grown  alloy  on  a  lattice  matched  substrate  (e.g.,  GalnAs  and  GaAsSb  grown  on 
InP,  and  GalnP  grown  on  GaAs),  the  energies  AE'  are  the  appropriate  excess  energies.  Since 
AE'  are  all  negative,  (see  Table  2)  these  alloys  will  be  in  the  ordered  states  at  a  temperature 
below  the  order-disorder  transition  temperature,  given  roughly  by  T0  =  I  AE'  I  /2k.  The  tem¬ 
peratures  T0  associated  with  the  tabulated  AE'  are  very  close  to  the  range  of  the  substrate  tem¬ 
peratures  at  which  the  LRO  were  found  to  exist  in  these  alloys.  This  may  be  a  mechanism  for 
the  LRO  in  these  alloys.  However,  this  mechanism  is  not  the  only  driving  force  for  the  LRO. 
For  the  GaAlAs  system,  the  calculated  excess  energy  is  essentially  zero,  because  GaAs  and 
AiAs  are  almost  lattice  matched.  Likewise,  we  are  also  puzzled  by  the  difficulty  in  the  growth 
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Table  2 


CALCULATED  VALUES  FOR  THREE  ORDERED  TERNARY  ALLOYS 


InGaAs 

InGaP 

GaSbAs 

CuAu  I 

cp 

CuAu  1 

cp 

CuAu  I 

cp 

d(A) 

2.537 

2.536 

2.454 

2.452 

2.550 

2.552 

d,  (A) 

2.602 

2.609 

2.523 

2.530 

2.627 

2.616 

d2  (A) 

2.476 

2.467 

2.389 

2.380 

2.447 

2.492 

Eb  (eV) 

-1.5762 

-1.5775 

-1.7431 

-1.7456 

-1.5451 

-1.5412 

AE  (meV) 

27 

22 

33 

23 

19 

35 

AE'  (meV) 

-92 

-87 

-106 

-96 

-130 

-120 

B  (101  'dynes/cm2) 

6.567 

6.599 

7.992 

8.028 

6.195 

6.111 

of  HgCdTe  and  the  experimental  findings  regarding  the  large  clustering  effects  despite  all  the 
mechanisms  that  we  know  of  that  give  a  very  small  mixing  enthalpy  and  indicate  that  HgCdTe 
should  behave  like  a  random  alloy. 

C.  Coulomb  Energies 

A  scheme  has  been  developed  to  calculate  the  configuration  dependence  of  the  Coulomb 
electron-electron  interactions  in  tetrahedrally  coordinated  pseudobinary  alloys  of  the  form 
A0  5B0.5C.9’10  Coulomb  energies  have  been  calculated  for  three  configurations:  a  random  alloy 
configuration,  an  ordered  sublattice  in  the  Cu-Au  structure,  and  a  spinodally  decomposed 
configuration.  The  Coulomb  energies,  driven  by  a  bond  polarity  mismatch,  must  be  considered 
as  a  contribution  to  the  mixing  enthalpy  in  semiconductor  alloys  in  addition  to  the  strain 
(driven  by  bond  length  mismatch)  and  chemical  terms  usually  considered.  The  results  of  the 
calculations  show  that  the  ordered  structures  have  a  lower  Coulomb  energy  than  the  random 
alloys.  This  is  because,  in  the  absence  of  long-range  order,  the  random  alloy  does  not  benefit 
from  a  Madelung  enhancement  of  the  electrostatic  energy  gained  from  a  coherent  sum  of  alter¬ 
nating  charges.  The  Coulomb  energy  terms  are  most  important  in  the  mixing  enthalpy  for  lat¬ 
tice  matched  systems  (e.g.,  AlGaAs,  where  experimental  evidence  of  ordering  has  been 
observed).11  Note  that  a  factor  of  two  errors  was  made  in  van  Schilfgaarde’s  papers9,10  and 
these  published  mixing  energies  are  a  factor  of  two  too  large. 
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Ill  VACANCIES  AND  SURFACE  SEGREGATION  IN  SEMICONDUCTORS 


A.  Vacancy  Formation  Energies 

Vacancies  can  affect  the  structural  properties  of  semiconductors  in  several  ways.  First, 
the  presence  of  vacancies  during  high-temperature  bulk  growth  may  form  vacancy  clusters  and 
precipitate  into  dislocation  loops.  Second,  the  presence  of  vacancies  is  important  in  both  sclf- 
and  impurity-diffusion  processes,  and  in  some  systems  are,  therefore,  related  to  the  efficiency 
of  annealing  and  doping.  Finally,  dislocation  motion  may  be  sensitive  to  the  presence  of 
vacancies,  via  climbing  processes,  whereby  vacancies  are  absorbed  or  emitted  and  through  an 
electrical  interaction  between  dislocations  and  vacancies,  as  in  the  case  of  charged  vacancies. 
Our  work  relating  to  the  calculation  of  the  total  energies  of  vacancy  formation  in  the  zinc- 
blende  compound  and  alloy  semiconductors  is  summarized  below;  general  conclusions  are 
included  in  Section  V. 

Vacancy  formation  energies  in  the  semiconductor  compounds  have  been  calculated  using  a 
tight-binding  cluster  model.12-14  Corrections  from  lattice  bond  length  relaxation  and  charge 
redistributions  about  the  vacancy  site  have  been  included.  Two  final  states  for  the  removed 
atom  have  been  considered  corresponding  to  a  free  atom  state  in  vacuum  and  to  the  removed 
atom  at  a  triply  bonding  site  on  a  (111)  surface,  and  a  doubly  bonding  site  on  a  (100)  surface. 
Cluster  convergence  has  been  examined  by  comparing  to  Green’s  function  calculation  that 
treats  interaction  to  infinite  range.  The  cluster  method  is  advocated  over  Green’s-function 
method,  because  it  allows  a  detailed  description  of  the  vacancy  formation  energies  in  alloys. 

Vacancy  formation  energies  in  pscudobinary  semiconductor  alloys  of  zinc-blende  structure 
have  been  calculated  using  a  modification  of  the  cluster  model  developed  for  the  com¬ 
pounds.15  Twelve  constituent  atoms  are  considered  (group  III:  Al,  Ga,  In;  group  V:  P,  As,  Sb; 
group  II:  Zn,  Cd,  Hg;  and  group  VI:  S,  Se,  Te),  and  cation  and  anion  extraction  energies  in  all 
III- V  and  II- VI  pseudobinary  alloys  formed  from  these  constituents  are  calculated.  The  extrac¬ 
tion  energies  are  found  to  be  sensitive  to  the  class  of  cluster  from  which  the  atom  is  removed, 
with  the  removal  of  a  cation  being  the  most  sensitive  to  the  alloy  environment  in  an  anion  sub¬ 
stituted  alloy,  and  likewise,  the  removal  of  an  anion  being  the  most  sensitive  to  the  alloy 
environment  in  a  cation  substituted  alloy.  In  most  random  alloys,  the  extraction  energies  are 
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found  to  vary  nonlinearly  with  alloy  composition,  and  most  often  bow  downwards,  in  some 
cases  reaching  a  minimum  at  a  finite  alloy  composition.  Such  nonlinearity  in  the  extraction 
energies  will  be  important  (e.g.,  in  the  interpretation  of  diffusion  and  doping  incorporation  in 
some  alloy  systems).  The  present  alloy  results  can  be  readily  extended  to  vapor  phase  finals 
states,  other  than  the  free  atom  assumed  for  the  extraction  energies  calculated  here,  by  includ¬ 
ing  the  appropriate  molecular  phase  energies  and  pressures.  Moreover,  the  extraction  energies 
can  be  directly  related  to  the  chemical  potentials  that  drive  interdiffusion  between  solids.  This 
is  especially  important  in  understanding  the  stability  of  epitaxial  layers  on  different  substrates 
and  of  superlattices. 

B.  Surface  Segregation 

The  driving  forces  for  surface  segregation  in  the  alloys,  whereby  the  surface  may  be  rich 
in  one  of  the  constituent  species,  are  related  to  the  vacancy  formation  energies  to  the  particular 
surface  of  interest.  We  have  studied  surface  segregation  in  several  systems,  and  a  summary  of 
this  work  is  given  below.  The  conclusions  that  we  have  drawn  are  discussed  in  Section  V. 

We  have  developed  a  theory  of  surface  segregation  for  the  pseudobinary  semiconductor 
alloys.16  Pseudobinary  alloys  of  the  form  Aj_xBxC  almost  always  have  different  surface  con¬ 
centrations  from  the  bulk  in  order  to  maintain  a  constant  chemical  potential  for  each  layer  in 
the  alloy.  The  enthalpy  responsible  for  segregation  is  the  difference  in  the  energies  for  moving 
an  A  or  B  atom  from  the  bulk  alloy  to  the  surface.  There  are  two  major  contributions  to  this 
energy:  (a)  a  bond-breaking  mechanism,  whereby  the  element  with  the  lowest  surface  energy 
segregates  to  the  top  and  (b)  strain  release,  where  the  dilute  element  in  the  compound  segre¬ 
gates  to  the  surface  to  alleviate  the  strain  energy  due  to  mismatch  of  the  AC  and  BC  bond 
lengths.  In  our  segregation  model,  the  free  energy  of  each  layer  is  calculated  in  the  regular 
and  quasichemical  approximatioas.  By  equating  the  chemical  potentials  of  each  successive 
layer  to  the  bulk,  the  composition  of  each  layer  is  obtained. 

The  degree  of  surface  segregation  for  the  alloys  Hgj_xCdxTe  (lattice  matched;  and 
Hg,_xZnxTc  (lattice  mismatched)  have  been  calculated.  Our  results  show  that  there  is  strong 
surface  enrichment  of  Hg  in  Hg)_xCdxTe,  while  there  is  less  surface  segregation  of  Hg  in 
Hgi-*ZnxTe  at  the  low  x  values  appropriate  for  infrared  applications.  Mercury  segregation  to 
the  surface  will  lower  the  band-gap  energy  or  may  turn  the  surface  into  a  semimetal,  thereby 
affecting  the  passivation  of  the  surface. 
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IV  MECHANICAL  PROPERTIES 
OF  SEMICONDUCTOR  COMPOUNDS  AND  ALLOYS 


A.  Elastic  Constants 

It  is  customary  to  use  a  two-parameter  force  constant  model  such  as  the  Keating  model17 
in  many  calculations  of  the  structure-related  properties  of  the  zinc-blende  material;  yet,  three 
independent  force  constants  are  necessary  to  describe  a  general  distortion  for  this  crystal  struc¬ 
ture.  We  have  developed  an  approximate,  but  ab  initio ,  theory  of  the  elastic  constants  in  the 
tetrahedral  semiconductors  that  provides  a  basis  for  the  generation  of  empirical  two-parameter 
force  constant  models18,1Q  A  relationship  between  the  three  independent  force  constants  has 
been  derived  that  is  independent  of  the  detailed  form  of  the  tight-binding  matrix  elements. 

We  have  developed  a  theory  for  the  internal  strain  energy  and  bulk  modulus  in  semi¬ 
conductor  alloys.20  The  internal  strain  energy  is  shown  to  be  governed  by  the  square  fluctua¬ 
tion  of  the  bond  lengths  and  an  effective  force  constant  that  couples  the  atomic  clusters  to  the 
environment.  Analytical  models  show  that  the  bulk  moduli  and  bond  lengths,  as  a  function  of 
concentration  of  these  alloys,  bow  only  slightly  below  the  straight-line  average.  This  work 
provides  a  basis  of  an  effective-medium  theory  for  calculating  the  input  cluster  energies  in  a 
statistical  theory. 

B.  Hardness 

Hardness  has  proven  to  be  a  useful  probe  of  the  mechanical  properties  of  the  brittle  semi¬ 
conductors.  Here  we  will  use  the  term  hardness  to  refer  specifically  to  Vickers’  hardness, 
unless  other  wise  noted.  In  the  Vickers’  hardness-measurement,  a  square  pyramidal  indenter  is 
used  and  the  hardness  number  is  given  by  the  applied  load  divided  by  the  area  of  the  indenta¬ 
tion  (i.e.,  units  of  pressure).  Hardness  has  been  found  to  be  an  intrinsic  property  of  the 
material,  because  it  is  relatively  independent  of  the  applied  load.  One  advantage  of  hardness 
measurements  for  semiconductors  is  that,  unlike  bending  tests,  only  small  samples  are  neces¬ 
sary  for  conventional  Vickers’  hardness-measurements,  or  for  nanoindenter  measurements,  rela¬ 
tively  thin  epitaxial  films  can  be  probed.  Additionally,  unlike  conventional  tests  used  to 
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measure  yield  stress,  hardness  measurements  can  be  made  at  room  temperature,  which  is  far 
below  the  usual  plastic  regime  for  most  semiconductors.  As  such,  the  hardness  measurement 
provides  a  convenient  and  usable  probe. 

The  question  remains,  though,  as  to  the  interpretation  of  the  hardness  measurement  in 
semiconductors;  that  is,  just  what  property  (or  properties)  of  a  semiconductor  are  we  measuring 
when  we  measure  hardness.  In  metals,  an  empirical  relationship  is  found  between  the  hard¬ 
ness,  H,  and  yield  stress,  Y : 

H  =  3  Y  . 

In  metals,  this  relationship  can  be  justified  of  the  basis  on  continuum  theory,  as  discussed  in 
McClintock  and  Argon.21  In  semiconductors,  such  a  simple  relationship  between  H  and  Y  is 
not  necessarily  appropriate  for  several  reasons.  During  deformation  in  metals,  many  slip 
planes  can  be  active  because  the  Peierls  barriers  for  dislocation  motion  in  most  directions  are 
low.  In  contrast,  because  the  bonds  in  semiconductors  are  strongly  covalent,  the  Peierls  bar¬ 
riers  are  high  and  dislocation  in  the  (111)  V$<TlO>  slip  system  dominate. 

To  date,  there  is  no  complete  quantitative  theory  of  hardness  in  the  semiconductors  in 
which  the  temperature  dependence,  photoplastic  effect,  and  the  alloy  hardening  effect  are 
included.  Sher  et  al.,22  proposed  a  model  of  hardness  for  the  semiconductor  compound  that 
gives  good  quantitative  agreement  with  experiment,  but  this  model  does  not  provide  an  expla¬ 
nation  for  several  of  the  observed  dependences  of  hardness.  The  Sher  model  of  hardness  in 
semiconductors  differs  from  more  conventional  interpretations  of  hardness;  In  the  Sher  model, 
the  hardness  is  dominated  by  dislocation-dislocation  interactions,  as  opposed  to  dislocation 
activation  and  motion  terms.  We  discuss  our  results  of  an  improved  quantitative  model  of 
hardness  below. 

I.  Slip  Systems 

In  the  Sher  hardness  model,  Vickers’  hardness  is  found  to  be  dominated  by  the  interaction 
energy  of  an  idealized  array  of  dislocations  that  have  been  generated  by  the  indenter.  In  their 
idealized  array,  no  account  was  taken  of  the  true  slip  systems  active  in  the  semiconductors. 
Experiments23  have  demonstrated  that,  in  Vickers’  hardness,  slip  occurs  primarily  on  the  {111} 
Vi  <U0>  glide  set,  where  the  three-fold  symmetry  of  slip  and  rosette  lines  occur  at  the  inter¬ 
section  of  the  {111}  planes  with  the  (111)  surface. 
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For  indentation  on  the  (1 1 1)  plane,  dislocations  can  glide  on  the  (1 1 1)  plane  parallel  to  the 
surface  or  on  one  of  the  three  other  {111}  planes  with  a  total  of  four  active  slip  planes. 
Although  the  detailed  analysis  differs  from  that  given  previously,22  the  contribution  to  the  hard¬ 
ness  from  the  interaction  energy  is  comparable  to  that  previously  calculated.  This  contribution 
to  H  is  directly  proportional  to  the  shear  coefficient. 

2.  Feierls  Energy 

The  Peierls  energy  is  difficult  to  calculate  precisely  because  of  dislocation  charge  effects 
and  reconstruction  at  the  dislocation  core.  In  the  context  of  the  hardness  measurement  we  wish 
calculate  the  Peierls  energy  to  access  the  importance  of  this  contribution  to  the  Vickers’  hard¬ 
ness  number.  Although  it  is  generally  agreed  that  dislocations  in  semiconductors  move 
through  the  generation  and  propagation  of  double  kinks,  in  the  hardness  measurement,  the 
region  about  the  indenter  is  grossly  plastically  deformed.  Because  the  dislocation  velocity  is 
low  at  room  temperature  (see  below),  the  large  dislocation  pile-up  model  proposed  by  Sher  et 
al.,22  may  be  appropriate.  If  the  dislocation  separation  is  small,  dislocation  motion  through 
kink  processes  will  be  suppressed,  and  the  dislocations  will  propagate  as  a  complete  unit. 

To  get  from  Configuration  A  to  Configuration  B  in  Figure  2,  we  must  break  a  row  of 
bonds.  Since  the  long-range  strain  fields  should  be  comparable  in  the  two  configurations  as 
well  as  for  intermediate  configurations,  the  Peierls  force  can  be  calculated  from  local  energy 
considerations  only.  The  energy  to  break  a  bond  at  the  dislocation  core  is  approximately  given 
by: 


Ub  =  2 VvJT  V32  +  2emet  -  V0  +  ^  [e£  -ehc] 


(1) 


where  V2  is  the  covalent  energy,  V3  the  ionic  energy,  £„,*  the  metallization  energy,  V0  the 
bond  overlap  energy,  and  and  the  hybrid  energy  for  the  anion  and  cation,  respectively. 
The  first  two  terms  in  Equation  (1)  account  for  the  loss  of  the  bonding  energy  of  the  two  elec¬ 
trons  in  the  breaking  bond,  the  third  term  accounts  for  regaining  the  repulsive  interaction 
energy  of  the  bond,  and  the  fourth  term  accounts  for  the  energy  gain  to  transfer  electrons  back 
from  the  cation  to  the  anion.  We  note  that  the  electron  orbitals  of  the  atoms  at  the  dislocation 
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core  are  left  in  the  sp3  hybrids  after  the  bond  breaking.  The  expression  in  Equation  (1) 
represents  a  theoretical  maximum  of  the  Peierls  energy,  since  no  reconstruction  at  the  core  has 
been  included. 

To  calculate  the  Peierls  energy  per  unit  length,  we  consider  a  primary  dislocation  in  the 
<1 10>  direction  in  a  zinc -blende  compound.  The  number  of  bonds  per  unit  length  in  <1 10>  is 
given  by  1/b  where  b  is  Burger’s  vector.  Thus,  the  Peierls  energy  per  unit  length  is  given  by 


(2) 


where  d  is  the  bond  length. 
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We  can  now  calculate  the  Peierls  force,  or  the  force  per  unit  length  necessary  to  move  a 
dislocation  over  the  potential  barrier,  as  illustrated  in  Figure  3.  The  Peieris  energy  is  related  to 
the  Peierls  force  through 


EP  =  FpL 


(3) 


where 


(4) 


is  the  distance  between  Configurations  A  and  Configuration  B.  Solving  for  FP  in  terms  of  Ub 
and  d,  we  arrive  at: 


Fp  = 


_3 

4 


Ub 

d2 


(5) 


or 


..JTli* 

V  2  8  d3 


(6) 


FIGURE  3  POTENTIAL  BARRIER  FOR  DISLOCATION  MOTION 
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Values  for  Ub  and  tP  are  summarized  in  Table  3. 


Table  3 


CALCULATED  PEIERLS  STRESS  AND  HARDNESS 
FOR  VARIOUS  ZINC  BLENDE  SEMICONDUCTORS. 
EXPERIMENTAL  HARDNESS  VALUES  SHOWN  FOR  COMPARISON. 


ub 

(ev) 

TP  2 
(kg/mm2) 

Hp  2 

(kg/mm  ) 

Hin. 

(kg/mm2) 

Hp  + 

(kg/mm2) 

Hexp 

(kg/mm2) 

c 

11.35 

23,300 

1940 

9244 

1 

Si 

5.95 

3430 

286 

1098 

Inf*  m 

1370 

Ge 

6.67 

3440 

286 

893 

1000 

Sn 

5.66 

1930 

161 

... 

mgm 

.... 

A1P 

6.04 

3410 

284 

___ 

... 

GaP 

6.15 

3500 

292 

903 

1195 

940 

InP 

5.78 

2640 

220 

548 

768 

520 

AlAs 

5.90 

3000 

250 

... 

_ 

505 

GaAs 

6.03 

3000 

256 

750 

1006 

580 

InAs 

5.64 

2340 

195 

469 

664 

430 

AlSb 

5.08 

2840 

237 

524 

761 

400 

GaSb 

5.36 

2180 

182 

553 

735 

450 

InSb 

4.92 

1670 

139 

365 

504 

230 

ZnS 

5.14 

3000 

250 

515 

765 

.... 

CdS 

4.77 

2750 

229 

288 

517 

.... 

HgS 

5.68 

2590 

216 

— 

— 

.... 

ZnSe 

4.92 

2500 

209 

462 

671 

137 

CdSe 

4.59 

1890 

157 

254 

411 

.... 

HgSe 

4.53 

1840 

154 

232 

386 

.... 

ZnTe 

4.42 

1720 

143 

374 

517 

82 

CdTe 

4.01 

1350 

113 

222 

335 

60 

HgTe 

3.99 

1360 

113 

230 

343 

25 
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Now  we  incorporate  the  Peierls  energy  into  the  hardness  model  for  low  temperature  where 
the  full  barrier  must  be  surmounted.  The  Sher  model  is  based  on  energy  considerations.  The 
Vickers’  hardness  number  is  given  by  the  applied  force  divided  by  the  area  of  indentation. 
Multiplying  the  numerator  and  denominator  by  h,  the  depth  of  indentation  we  have 
H  =  E/(W2h),  where  E  =  F  h  is  the  energy  of  indentation,  and  h  is  the  depth  of  the  indentation. 
Including  the  interaction  energy  only  we  have: 


H  =  Hmt  = 


G  cot8 
6rc(l-v) 


-ln(^)  +  i  +  sin’-f 


where  0  is  one  half  the  indenter  angle.  To  include  the  Peierls  energy,  we  consider  the  total 
energy  necessary  to  move  the  dislocations  from  their  initial  to  final  positions  in  the  idealized 
model.  The  hardness  is  then  given  by: 


H  =  HmI  +  HP 


where 


HP 


EP 

W2h 


is  the  Peierls  contribution  to  the  hardness,  and  EP  is  the  total  Peierls  energy  expended.  The 
total  length  of  dislocation  to  be  moved  is  calculated  to  be: 


.  1  WJ 

Ur~6  ? 


cos20 


The  total  Peierls  energy  is  given  by: 


<J„.  1  Usw- 


COS2© 
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Thus,  we  have: 


HP 


1  Ub  q  •  Q 

- -cos0  sin0 

3  b3 


For  0  =  45°: 


Values  of  HP  arc  summarized  in  Table  3.  Several  features  of  HP  should  be  noted.  First, 
we  have  used  a  zero  temperature  value  of  the  Peierls  energy.  Because  hardness  measurements 
are  typically  done  at  room  temperature,  one  should  take  the  thermal  energy  into  account,  which 
will  reduce  the  values  of  HP  from  those  listed  in  Table  3.  Also  shown  in  Table  3  are 
Him.  Him+Hp,  (the  best  theoretical  estimate  for  H),  and  Hexp.  Note  that,  like  Hi,,,,  Hp  is 
independent  of  the  applied  load,  in  agreement  with  experiment.  Also  note  that  Hp  improves 
the  agreement  between  theory  and  experiment  for  the  hard,  nonpolar  materials.  For  the  softer, 
more  ionic  materials,  H  is  overestimated  by  the  theory.  We  have  begun  to  study  the  charge 
and  polarity  dependence  of  the  dislocation  motion,  and  expect  that  the  over  estimation  of  H  is 
because  of  neglect  of  dislocation  velocity  effects  and  their  temperature  dependence,  as  dis¬ 
cussed  below. 

3.  Temperature  Dependence 

We  have  begun  to  develop  a  theory  of  the  explicit  temperature  dependence  of  hardness. 
Here  we  summarize  the  experimental  results  and  discuss  a  tentative  theory  of  the  temperature 
dependence. 

Several  recent  studies  on  the  temperature  dependence  of  hardness  serve  to  illustrate  the 
behavior.  Results  for  GaAs  and  Ge  are  shown  in  Figure  4.  The  (111)  and  (100)  faces  of 
GaAs  have  been  examined  by  Hirsch  et  al.,23  and  Guruswamy  et  al.,24  respectively.  Results 
for  the  Knoop  hardness  on  the  (100)  face  of  n-type  Ge  are  also  shown.25  The  (100)  face  of 
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FIGURE  4  TEMPERATURE  DEPENDENCE  OF  HARDNESS  FOR  GaAs  AND  Ge 


GaAs  and  the  Ge  show  a  definite  temperature  dependence  with  a  relatively  temperature 
independent  region  for  T  <  450K  and  an  exponential  temperature  dependence  for  T  >  550  K: 


u 

H  =  Hc  e  kT 


with 


U  S  0.24  eV 

for  (100)  GaAs.  The  results  for  GaAs  (111)  appear  to  follow  similar  behavior. 

The  temperature  dependence  of  hardness  suggests  that  two  different  mechanisms  may 
determine  hardness  in  the  two  temperature  regimes.  At  low  temperature,  the  hardness  is  nearly 
independent  of  temperature  and  may  be  limited  primarily  by  dislocation  interactions.  Disloca¬ 
tion  mobility  is  low  at  low  temperatures,  and  the  tendency  for  dislocation  pile-up  is  high.  At 
elevated  temperature,  the  dislocation  mobility  is  increased,  so  that  dislocations  move  more 
readily  under  an  applied  stress.  Therefore,  at  higher  temperatures,  dislocation  pile-up  is 
reduced  and  the  hardness  is  limited  by  lattice  friction,  which  shows  a  strong  temperature 
dependence.  We  are  currently  investigating  a  quantitative  theory  of  the  temperature  depen¬ 
dence  of  hardness,  which  also  includes  the  dissimilar  velocities  of  the  a  and  p  dislocations. 


C.  Cleavage  and  Fracture 

One  of  the  basic  ingredients  to  the  calculation  of  the  fracture  properties  in  real  materials  is 
the  cleavage  energy.  For  ideal  brittle  fracture,  i.e.,  fracture  unaccompanied  by  crack  tip  blunt¬ 
ing,  dislocation  generation,  and  phonon  generation,  the  cleavage  energy  provides  a  good  esti¬ 
mate  of  the  fracture  energy. 

A  simple  model  is  usually  used  to  calculate  the  cleavage  energy,  e,  based  on  the  energy  to 
separate  two  halves  of  a  semiconductor  crystal.  That  is: 


a(x)dx 
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where  o(x)  is  the  restoring  force  per  unit  area  for  the  crystal  halves  separated  by  a  distance  x, 
where  x  =  0  is  the  equilibrium  position.  We  have  begun  to  develop  a  method  to  calculate 
o(x),  based  on  the  cleavage  of  a  crystal  of  infinite  extent  This  involved  the  calculation  of  the 
total  energy  of  two  semi-infinite  crystals  as  a  function  of  their  separation.  Based  on  a  tight- 
binding  theory,  the  cleavage  energy,  lime(x),  can  be  written  as: 


E(X)  =  Eel  +  V0  +  Ceouj 


where  ecl  is  the  change  in  the  electron  energies,  VD  is  the  sum  of  the  repulsive  pair  energies 
between  the  atom  pairs  across  the  interface,  and  is  the  Coulomb  energy,  important  for 
cleavage  on  polar  surfaces.  It  is  that  is  difficult  to  calculate  accurately;  can  be  calcu¬ 
lated  systematically  using  the  theory  of  the  interface  Green’s  function  we  have  developed.  A 
summary  of  this  method  follows. 

The  theory  of  surface  Green’s  function  has  been  developed  for  general  layered  structures, 
where,  for  the  present  case  of  cleavage,  the  layers  are  considered  parallel  to  the  cleavage  sur¬ 
face.  For  a  given  crystal  surface  direction,  here  assumed  to  be  the  {110}  cleavage  plane,  all 
the  Green’s  functions  are  ^/-dependent.  That  is: 

£ei  =  Z 

where 

£*!<£//)  =  J  ApOc/z.EiEdE  +  £  ANjCj 


where  E  is  the  energy  and  e*  is  the  Fermi  energy,  and  where 


Ap04.E)  =  I  Apn0c//tE) 

n 
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where  n  is  the  layer  index  and  varies  from  -=»  to  +°°,  and  the  e*  are  possible  localized  states  at 
the  surface.  The  cleavage  is  assumed  to  occur  between  the  layers  n  =  0  and  n  =  1.  Apn  (£„, E) 
is  given  by: 


Apn(?/;f)  =  -  J^ImGnn  -  Imgnnj  . 

Gnn  and  gn,,  are  the  interface  and  surface  Green’s  functions,  respectively.  These  Green’s  func¬ 
tions  are  both  %,  and  energy  dependent.  The  calculation  of  these  Green’s  functions  are  straight 
forward  and  work  has  begun  on  setting  up  the  problem  for  the  {110)  surfaces  in  zinc-blende 
semiconductor  compounds.  The  calculation  of  this  cleavage  energy  is  simplified  by  the  fact 
that  the  {110}  surfaces  are  nonpolar. 
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V  CONCLUSIONS 


The  major  conclusions  of  our  work  are  summarized  below: 

.  Semiconductor  alloys  are  never  truly  random,  and  deviations  from  randomness  have 
been  calculated  for  many  III-V  and  II-V  compounds. 

.  Deviations  from  random  alloys  are  driven  by  excess  energies  due  to  bond-length 
differences,  chemical  differences,  and  polarity-induced  Coulomb-energy  differences 
between  the  constituents,  and  by  externally  applied  coherent  strain  fields  and  tempera¬ 
ture  gradients.  The  correlation  status  might  be  growth-method  dependent. 

.  Coulomb  energies  for  ordered  structures  are  smaller  than  for  random  alloys  because  of 
a  Madelung  enhancement  in  the  ordered  structures. 

.  Bond  polarities  have  been  calculated  in  excellent  agreement  with  experiment 

.  Vacancy  formation  energies  in  the  semiconductor  compounds  are  found  to  be  surface- 
orientation  dependent. 

•  Vacancy  formation  energies  in  the  semiconductor  alloys  are  found  to  vary  nonlinearly 
with  the  near-neighbor  alloy  composition.  Several  possible  consequences  of  this  arc: 

-  Composition-  and  correlation-state  dependent  vacancy  chemical  potentials:  multiple- 
vacancy-activation  energies  and  entropies,  and  anomalous  concentrations. 

-  Anomalous  diffusion:  multiple-diffusion-activation  energies,  and  entropies;  parallel- 
path  diffusion  through  percolating  channels. 

.  A  new  effective  medium  theory  predicts  an  approximately  linear  variation  of  the  bulk 
modules  with  alloy  composition. 

.  A  foundation  for  the  reduction  from  three  to  two-parameter  force-constant  models 
intetrahedrally  coordinated  has  been  established  via  an  ab  initio  theory. 

.  The  model  of  hardness  has  been  improved  by  including  the  Peieris  energies  into  the 
model. 

.  A  theory  of  surface  segregation  has  been  initiated  that  shows  surfaces: 

-  Above  Tc  (the  critical  spinodal  temperature)  that  are  rich  in  one  component  and 
where  the  concentration  returns  to  the  bulk  concentration  in  a  few  atomic  layers. 

-  Below  Tc  in  which  preliminary  results  show  that  the  concentration  oscillates  back  to 
the  bulk  value  over  many  lattice  spacings.  This  is  a  superlattice. 


21 


REFERENCES 


1.  A.  Sher,  M.  van  Schilfgaarde,  A.-B.  Chen,  and  W.  Chen,  “Quasichemical  approximation 
in  binary  alloys,”  Phys.  Rev.  B,  Vol.  36,  p.  4279,  1987. 

2.  A.  Sher,  A.-B.  Chen,  M.  van  Schilfgaarde,  and  M.A.  Berding,  “Semiconductor  alloys: 
structural  property  engineering,”  Proc.  1st  International  SAMPE  Electronics  Conference 
(23-25  June  1987,  Sunnyvale,  California),  p.  323,  1987. 

3.  A.  Sher,  M.A.  Berding,  S.  Krishnamurthy,  M.  van  Schilfgaarde,  A.-B.  Chen,  and  W. 
Chen,  “Structure-property  relationships  in  semiconductor  alloys,”  Mat.  Res.  Soc.  Symp. 
Proc.  Vol.  90,  91  (1987c) 

4.  A.  Sher,  M.A.  Berding,  M.  van  Schilfgaarde,  A.-B.  Cnen,  and  R.S.  Patrick,  “Modeling  of 
mechanical  properties  of  II- VI  materials,"  J.  Cryst.  Growth,  Vol.  86,  p.  15,  1988. 

5.  R.S.  Patrick,  A.-B.  Chen,  A.  Sher,  and  M.A.  Berding,  “Phase  diagrams  and  microscopic 
structures  of  (Hg.Cd)Te,  (Hg,Zn)Te,  and  (Cd,Zn)Te  alloys,”  J.  Vac.  Sci.  Tech.  A,  Vol.  6, 
p.  2643,  1988. 

6.  M.A.  Berding,  A.  Sher,  A.-B.  Chen,  and  W.E.  Miller,  “Structural  properties  of  bismuth¬ 
bearing  semiconductor  alloys,”  J.  Appl.  Phys.,  Vol.  63,  p.  107,  1988. 

7.  R.S.  Patrick,  A.-B.  Chen,  and  A.  Sher,  “Phase  diagrams  and  local  correlations  in  pseudo- 
binary  alloys,”  Phys.  Rev.  B,  Vol.  36,  p.  6585,  1987. 

8.  G.B.  Stringfellow,  and  P.E.  Greene,  J.  Phys.  Chem.  Solids,  Vol.  30,  p.  1779,  1969. 

9.  M.  vanSchilfgaarde,  A.-Ban  Chen,  and  A.  Sher,  "Coulomb  energy  in  pseudobinary 
alloys,”  Phys.  Rev.  Lett.,  Vol.  57,  p.  1149,  1986. 

10.  M.  van  Schilfgaarde,  A.-Ban  Chen,  and  A.  Sher,  “Configuration  dependence  of  coulomb 
interactions  in  pseudobinary  alloys,”  Proceedings  of  the  1986  7th  International  Confer¬ 
ence  on  Ternary  and  Multinary  Compounds,  10-12  September  Snowmass,  Colorado,  p. 
411,  1987. 


22 


11.  T.S.  Kuan,  T.F.  Kuech,  W.I.  Wang,  and  E.L.  Wilkie,  Phys.  Rev.  Lett.,  Vol.  54,  p.  201, 
1985. 

12.  M.A.  Bending,  A.  Sher,  A.-B.  Chen,  “Vacancy  formation  energies  in  II-VI  semi¬ 
conductors,”  J.  Vac.  Sci.  Technol.  A,  Vol.  5,  p.  3009,  1987. 

13.  M.A.  Berding,  A.-B.  Chen,  and  A.  Sher,  “HgCdTe  versus  HgZnTe:  Electronic  properties 
and  vacancy  formation  energies,”  Mat.  Res.  Soc.  Symp.  Proc.,  Vol.  90,  p.  127,  1987. 

14.  M.A.  Berding,  A.  Sher,  A.-B.  Chen,  “Vacancy  formation  energy  in  semiconductor  com¬ 
pounds,”  submitted  to  Physical  Review  B. 

15.  M.A.  Berding,  A.  Sher,  A.-B.  Chen,  “Vacancy  formation  energy  in  semiconductor 
alloys,”  in  preparation. 

16.  R.S.  Patrick,  A.-B.  Chen,  A.  Sher,  and  M.A.  Berding,  “Surface  Segregation  in  Pseudo¬ 
binary  Alloys,”  in  preparation. 

17.  P.N.  Keating,  “Effect  of  invariance  requirements  on  the  elastic  strain  energy  of  crystals 
with  application  to  the  diamond  structure,”  Phys.  Rev.,  Vol.  145,  p.  637,  1966. 

18.  M.  van  Schilfgaarde,  and  A.  Sher,  “Tight-binding  theory  and  elastic  constants,”  Phys. 
Rev.  B,  Vol.  36,  p.  4375,  1987. 

19.  M.  van  Schilfgaarde,  and  A.  Sher,  “Tight-binding  theory  of  force  constant  models,”  Appl. 
Phys.  Lett,  Vol.  51,  p.  175,  1987. 

20.  A.-B.  Chen,  A.  Sher,  and  M.A.  Berding,  “Semiconductor  alloy  theory:  Internal  strain 
energy  and  bulk  modulus,”  Phys.  Rev.  B,  Vol.  37,  p.  6285,  1988. 

21.  F.A.  McClintock,  and  A.S.  Argon,  Mechanical  Behavior  of  Materials,  Addison-Wesley, 
1966. 

22.  A.  Sher,  A.-B.  Chen,  and  W.E.  Spicer,  “Dislocation  energies  and  hardness  of  semi¬ 
conductors,”  Appl.  Phys.  Lett.,  Vol.  46,  p.  54,  1985. 

23.  P.B.  Hirsch,  P.  Pirouz,  S.G.  Roberts,  and  P.D.  Warren,  “Indentation  Plasticity  and  Polar¬ 
ity  of  Hardness  on  (111)  Faces  of  GaAs,”  Philos.  Mag.  B,  Vol.  52,  p.  759,  1985. 


23 


24.  S.  Guruswamy,  J.P.  Hirth,  K.T.  Faber,  “High  Temperature  Hardness  of  Ga,_xInxAs,”  J 
Appl.  Phys.,  Vol.  60,  p.  4136,  1986. 

25.  S.G.  Roberts,  P.D.  Warren,  and  P.B.  Hirsch,  “Knoop  Hardness  Anisotropy  on  (001) 
Faces  of  Germanium  and  GaAs,”  J.  Mater.  Res.,  Vol.  1,  p.  162,  1986. 


24 


BIBLIOGRAPHY 


M.  vanSchilfgaarde,  A.-Ban  Chen,  and  A.  Sher,  “Coulomb  energy  in  pseudobinary  alloys,” 
Phys.  Rev.  Lett.,  Vol.  57,  p.  1149,  1986. 

M.A.  Berding,  A.  Sher,  A.-B.  Chen,  “Vacancy  formation  energies  in  II-VI  semiconductors,”  J. 
Vac.  Sci.  Technol.  A,  Vol.  5,  p.  3009,  1987. 

M.A.  Berding,  A.  Sher,  and  A.-B.  Chen,  “Polarity  in  Semiconductor  Compounds,”  Phys.  Rev. 
B,  Vol.  36,  p.  7433,  1987. 

M.A.  Berding,  A.-B.  Chen,  and  A.  Sher,  “HgCdTe  versus  HgZnTc:  Electronic  properties  and 
vacancy  formation  energies,”  Mat.  Res.  Soc.  Symp.  Proc.,  Vol.  90,  p.  127,  1987. 

A.-B.  Chen,  M.  van  Schilfgaarde,  S.  Krishnamurthy,  M.A.  Berding,  A.  Sher,  “Alloy  electronic 
structure  and  statistics,”  Proceedings  of  the  1986  7th  International  Conference  on  Ternary 
and  Multinary  Compounds,  10-12  September  Snowmass,  Colorado,  p.  377,  1987. 

R.S.  Patrick,  A.-B.  Chen,  and  A.  Sher,  "Phase  diagrams  and  local  correlations  in  pseudobinary 
alloys,”  Phys.  Rev.  B,  Vol.  36,  p.  6585,  1987. 

A.  Sher,  M.A.  Berding,  S.  Krishnamurthy,  M.  van  Schilfgaarde,  A.-B.  Chen,  and  W.  Chen, 
“Structure-property  relationships  in  semiconductor  alloys,”  Mat.  Res.  Soc.  Symp.  Proc., 
Vol.  90.  p.  91,  1987. 

A.  Sher,  A.-B.  Chen,  M.  van  Schilfgaarde,  and  M.A.  Berding,  “Semiconductor  alloys:  struc¬ 
tural  property  engineering,"  Proc.  1st  International  SAMPLE  Electronics  Conference.  Sun¬ 
nyvale,  California,  p.  323,  1987. 

A.  Sher,  M.  van  Schilfgaarde,  A.-B.  Chen,  and  W.  Chen,  “Quasichemical  approximation  in 
binary  alloys,”  Phys.  Rev.  B,  Vol.  36,  p.  4279,  1987. 


25 


M.  van  Schilfgaarde,  A.-Ban  Chen,  and  A.  Sher,  “Configuration  dependence  of  coulomb 
interactions  in  pseudobinary  alloys”  Proceedings  of  the  1986  7th  International  Conference 
on  Ternary  and  Multinary  Compounds,  10-12  September  Snowmass,  Colorado,  p.  411, 
1987. 

M.  van  Schilfgaarde  and  A.  Sher,  “Tight-binding  theory  and  elastic  constants,”  Phys.  Rev.  B, 
Vol.  36,  p.  4375,  1987. 

M.  van  Schilfgaarde  and  A.  Sher,  “Tight-binding  theory  of  force  constant  models,”  Appl. 
Phys.  Lett.,  Vol.  51,  p.  175,  1987. 

M.A.  Berding,  A.  Sher,  A.-B.  Chen,  and  W.E.  Miller,  “Structural  properties  of  bismuth¬ 
bearing  semiconductor  alloys, q  J.  Appl.  Phys.,  Vol.  63,  p.  107,  1988. 

A.-B.  Chen,  A.  Sher,  and  M.A.  Berding,  “Semiconductor  alloy  theory:  Internal  strain  energy 
and  bulk  modulus,”  Phys.  Rev.  B,  Vol.  37,  p.  6285,  1988. 

R.S.  Patrick,  A.-B.  Chen,  A.  Sher,  and  M.A.  Berding,  “Phase  diagrams  and  microscopic  struc¬ 
tures  of  (Hg,Cd)Te,  (Hg,Zn)Te.  and  (Cd,Zn)Te  alloys,”  /.  Vac.  Sci.  Tech.  A,  Vol.  6,  p. 
2643,  1988.  ' 

A.  Sher,  M.A.  Berding,  A.-B.  Chen,  and  R.  Patrick,  “Stability  of  Alloy  Semiconductors:  Com¬ 
mon  Tellurium  Anion  II-IV  Semiconductor,”  presented  at  “The  Future  of  Small  Gap  II- 
VI  Semiconductors,”  Univ.  de  Liege,  Belgium,  September,  1988. 

A.  Sher,  M.A.  Berding,  M.  van  Schilfgaarde,  A.-B.  Chen,  and  R.S.  Patrick,  “Modeling  of 
mechanical  properties  of  II-VI  materials,”  /.  Cryst.  Growth,  Vol.  86,  p.  15,  1988. 

M.A.  Berding,  A.  Sher,  A.-B.  Chen,  “Vacancy  formation  energy  in  semiconductor  com¬ 
pounds,”  submitted  to  Physical  Review  B. 

A.-B.  Chen,  Y.-M.  Lai-Hsu,  and  W.  Chen,  “Difference-equation  approach  to  the  electronic 
structures  of  surfaces,  interfaces  and  superlattices,”  submitted  to  Physical  Review  B. 

M.A.  Berding,  A.  Sher,  A.-B.  Chen,  "Vacancy  formation  energy  in  semiconductor  alloys,”  in 
preparation. 


26 


R.S.  Patrick,  A.-B.  Chen,  A.  Sher,  and  M.A.  Berding,  “Surface  Segregation  in  pseudobinary 
alloys,”  in  preparation. 


27 


